The evolving SARS-CoV-2 epidemic in Africa: Insights from rapidly expanding genomic surveillance. (preprint)

Investment in Africa over the past year with regards to SARS-CoV-2 genotyping has led to a massive increase in the number of sequences, exceeding 100,000 genomes generated to track the pandemic on the continent. Our results show an increase in the number of African countries able to sequence within their own borders, coupled with a decrease in sequencing turnaround time. Findings from this genomic surveillance underscores the heterogeneous nature of the pandemic but we observe repeated dissemination of SARS-CoV-2 variants within the continent. Sustained investment for genomic surveillance in Africa is needed as the virus continues to evolve, particularly in the low vaccination landscape. These investments are very crucial for preparedness and response for future pathogen outbreaks.

A year of genomic surveillance reveals how the SARS-CoV-2 pandemic unfolded in Africa (preprint)

The progression of the SARS-CoV-2 pandemic in Africa has so far been heterogeneous and the full impact is not yet well understood. Here, we describe the genomic epidemiology using a dataset of 8746 genomes from 33 African countries and two overseas territories. We show that the epidemics in most countries were initiated by importations, predominantly from Europe, which diminished following the early introduction of international travel restrictions. As the pandemic progressed, ongoing transmission in many countries and increasing mobility led to the emergence and spread within the continent of many variants of concern and interest, such as B.1.351, B.1.525, A.23.1 and C.1.1. Although distorted by low sampling numbers and blind-spots, the findings highlight that Africa must not be left behind in the global pandemic response, otherwise it could become a breeding ground for new variants.

Genomic epidemiology of the SARS-CoV-2 epidemic in Zimbabwe: Role of international travel and regional migration in spread (preprint)

Zimbabwe reported its first case of SARS-Cov-2 infection in March 2020, and case numbers increased to more than 8,099 to 16th October 2020. An understanding of the SARS-Cov-2 outbreak in Zimbabwe will assist in the implementation of effective public health interventions to control transmission. Nasopharyngeal samples from 92,299 suspected and confirmed COVID-19 cases reported in Zimbabwe between 20 March and 16 October 2020 were obtained. Available demographic data associated with those cases identified as positive (8,099) were analysed to describe the national breakdown of positive cases over time in more detail (geographical location, sex, age and travel history). The whole genome sequence (WGS) of one hundred SARS-CoV-2-positive samples from the first 120 days of the epidemic in Zimbabwe was determined to identify their relationship to one another and WGS from global samples. Overall, a greater proportion of infections were in males (55.5%) than females (44.85%), although in older age groups more females were affected than males. Most COVID-19 cases (57 %) were in the 20-40 age group. Eight lineages, from at least 25 separate introductions into the region were found using comparative genomics. Of these, 95% had the D614G mutation on the spike protein which was associated with higher transmissibility than the ancestral strain. Early introductions and spread of SARS-CoV-2 were predominantly associated with genomes common in Europe and the United States of America (USA), and few common in Asia at this time. As the pandemic evolved, travel-associated cases from South Africa and other neighbouring countries were also recorded. Transmission within quarantine centres occurred when travelling nationals returning to Zimbabwe. International and regional migration followed by local transmission were identified as accounting for the development of the SARS-CoV-2 epidemic in Zimbabwe. Based on this, rapid implementation of public health interventions are critical to reduce local transmission of SARS-CoV-2. Impact of the predominant G614 strain on severity of symptoms in COVID-19 cases needs further investigation.

Protein scaffold-based multimerization of soluble ACE2 efficiently blocks SARS-CoV-2 infection in vitro (preprint)

Soluble ACE2 (sACE2) decoy receptors are promising agents to inhibit SARS-CoV-2 as they are not affected by common escape mutations in viral proteins. However, their success may be limited by their relatively poor potency. To address these challenges, we developed a highly active multimeric sACE2 decoy receptor via a SunTag system that could neutralize both pseudoviruses bearing SARS-CoV-2 spike protein and SARS-CoV-2 clinical isolates. This fusion protein demonstrated a neutralization efficiency nearly 250-fold greater than monomeric sACE2. SunTag in combination with a more potent version of sACE2 achieved near complete neutralization at a sub-nanomolar range, which is comparable with clinical monoclonal antibodies. We demonstrate that this activity is due to greater occupancy of the multimeric decoy receptors on Spike protein as compared to monomeric sACE2. Overall, these highly potent multimeric sACE2 decoy receptors offer a promising treatment approach against SARS-CoV-2 infections including its novel variants.

Chitinase 3-like-1 is a Therapeutic Target That Mediates the Effects of Aging in COVID-19 (preprint)

COVID-19 is caused by the SARS-CoV-2 (SC2) virus and is more prevalent and severe in the elderly and patients with comorbid diseases (CM). Because chitinase 3-like-1 (CHI3L1) is induced during aging and CM, the relationships between CHI3L1 and SC2 were investigated. Here we demonstrate that CHI3L1 is a potent stimulator of the SC2 receptor ACE2 and viral spike protein priming proteases (SPP), that ACE2 and SPP are induced during aging and that anti-CHI3L1, kasugamycin and inhibitors of phosphorylation, abrogate these ACE2- and SPP- inductive events. Human studies also demonstrated that the levels of circulating CHI3L1 are increased in the elderly and patients with CM where they correlate with COVID-19 severity. These studies demonstrate that CHI3L1 is a potent stimulator of ACE2 and SPP; that this induction is a major mechanism contributing to the effects of aging during SC2 infection and that CHI3L1 coopts the CHI3L1 axis to augment SC2 infection. CHI3L1 plays a critical role in the pathogenesis of and is an attractive therapeutic target in COVID-19.

SARS-CoV-2 spike downregulates tetherin to enhance viral spread (preprint)

The antiviral restriction factor, tetherin, blocks the release of several different families of enveloped viruses, including the Coronaviridae. Tetherin is an interferon-induced protein that forms parallel homodimers between the host cell and viral particles, linking viruses to the surface of infected cells and inhibiting their release. We demonstrate that SARS-CoV-2 downregulates tetherin to aid its release from cells, and investigate potential proteins involved in this process. Loss of tetherin from cells caused an increase in SARS-CoV-2 viral titre. We find SARS-CoV-2 spike protein to be responsible for tetherin downregulation, rather than ORF7a as previously described for the 2002-2003 SARS-CoV. We instead find ORF7a to be responsible for Golgi fragmentation, and expression of ORF7a in cells recapitulates Golgi fragmentation observed in SARS-CoV-2 infected cells.

Distinct mutations and lineages of SARS-CoV-2 virus in the early phase of COVID-19 global pandemic and subsequent global expansion (preprint)

A novel coronavirus, SARS-CoV-2, has caused over 8538 million cases and over 1.8 1 million deaths worldwide since it occurred twelve months ago in Wuhan, China. Here we first analyzed 4,013 full-length SARS-CoV-2 genomes from different continents over a 14-week timespan since the outbreak in Wuhan, China. 2,954 unique nucleotide substitutions were identified with 31 of the 4,013 genomes remaining as ancestral type, and 952 (32.2%) mutations recurred in more than one genome. A viral genotype from the Seafood Market in Wuhan featured with two concurrent mutations was the dominant genotype (80.9%) of the pandemic. We also identified unique genotypic compositions from different geographic locations, and time-series viral genotypic dynamics in the early phase that reveal transmission routes and subsequent expansion. We also used the same approach to analyze 261,350 full-length SARS-CoV-2 genomes from the world over 12 months since the outbreak (i.e. all the available viral genomes in the GISAID database as of 25 December 2020). Our study indicates the viral genotypes can be utilized as molecular barcodes in combination with epidemiologic data to monitor the spreading routes of the pandemic and evaluate the effectiveness of control measures.

COVID-19 virtual patient cohort reveals immune mechanisms driving disease outcomes (preprint)

To understand the diversity of immune responses to SARS-CoV-2 and distinguish features that predispose individuals to severe COVID-19, we developed a mechanistic, within-host mathematical model and virtual patient cohort. Our results indicate that virtual patients with low production rates of infected cell derived IFN subsequently experienced highly inflammatory disease phenotypes, compared to those with early and robust IFN responses. In these in silico patients, the maximum concentration of IL-6 was also a major predictor of CD8+ T cell depletion. Our analyses predicted that individuals with severe COVID-19 also have accelerated monocyte-to-macrophage differentiation that was mediated by increased IL-6 and reduced type I IFN signalling. Together, these findings identify biomarkers driving the development of severe COVID-19 and support early interventions aimed at reducing inflammation.

RNA-protein interaction analysis of SARS-CoV-2 5'- and 3'-untranslated regions identifies an antiviral role of lysosome-associated membrane protein-2 (preprint)

Severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) is a positive-strand RNA virus. Viral genome is capped at the 5'-end, followed by an untranslated region (UTR). There is poly-A tail at 3'-end, preceded by an UTR. Self-interaction between the RNA regulatory elements present within 5'- and 3'-UTRs as well as their interaction with host/virus-encoded proteins mediate the function of 5'- and 3'-UTRs. Using RNA-protein interaction detection (RaPID) assay coupled to liquid chromatography with tandem mass-spectrometry, we identified host interaction partners of SARS-CoV-2 5'- and 3'-UTRs and generated an RNA-protein interaction network. By combining these data with the previously known protein-protein interaction data proposed to be involved in virus replication, we generated the RNA-protein-protein interaction (RPPI) network, likely to be essential for controlling SARS-CoV-2 replication. Notably, bioinformatics analysis of the RPPI network revealed the enrichment of factors involved in translation initiation and RNA metabolism. Lysosome-associated membrane protein-2a (Lamp2a) was one of the host proteins that interact with the 5'-UTR. Further studies showed that Lamp2 level is upregulated in SARS-CoV-2 infected cells and overexpression of Lamp2a and Lamp2b variants reduced viral RNA level in infected cells and vice versa. In summary, our study provides an useful resource of SARS-CoV-2 5'- and 3'-UTR binding proteins and reveal the antiviral function of host Lamp2 protein.

Human cell-dependent, directional, time-dependent changes in the mono- and oligonucleotide compositions of SARS-CoV-2 genomes (preprint)

Background: When a virus that has grown in a nonhuman host starts an epidemic in the human population, human cells may not provide growth conditions ideal for the virus. Therefore, the invasion of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), which is usually prevalent in the bat population, into the human population is thought to have necessitated changes in the viral genome for efficient growth in the new environment. In the present study, to understand host-dependent changes in coronavirus genomes, we focused on the mono- and oligonucleotide compositions of SARS-CoV-2 genomes and investigated how these compositions changed time-dependently in the human cellular environment. We also compared the oligonucleotide compositions of SARS-CoV-2 and other coronaviruses prevalent in humans or bats to investigate the causes of changes in the host environment. Results: Time-series analyses of changes in the nucleotide compositions of SARS-CoV-2 genomes revealed a group of mono- and oligonucleotides whose compositions changed in a common direction for all clades, even though viruses belonging to different clades should evolve independently. Interestingly, the compositions of these oligonucleotides changed towards those of coronaviruses that have been prevalent in humans for a long period and away from those of bat coronaviruses. Conclusions: Clade-independent, time-dependent changes are thought to have biological significance and should relate to viral adaptation to a new host environment, providing important clues for understanding viral host adaptation mechanisms.

Large scale sequencing of SARS-CoV-2 genomes from one region allows detailed epidemiology and enables local outbreak management (preprint)

The COVID-19 pandemic has spread rapidly throughout the world. In the UK, the initial peak was in April 2020; in the county of Norfolk (UK) and surrounding areas, which has a stable, low-density population, over 3,200 cases were reported between March and August 2020. As part of the activities of the national COVID-19 Genomics Consortium (COG-UK) we undertook whole genome sequencing of the SARS-CoV-2 genomes present in positive clinical samples from the Norfolk region. These samples were collected by four major hospitals, multiple minor hospitals, care facilities and community organisations within Norfolk and surrounding areas. We combined clinical metadata with the sequencing data from regional SARS-CoV-2 genomes to understand the origins, genetic variation, transmission and expansion (spread) of the virus within the region and provide context nationally. Data were fed back into the national effort for pandemic management, whilst simultaneously being used to assist local outbreak analyses. Overall, 1,565 positive samples (172 per 100,000 population) from 1,376 cases were evaluated; for 140 cases between two and six samples were available providing longitudinal data. This represented 42.6% of all positive samples identified by hospital testing in the region and encompassed those with clinical need, and health and care workers and their families. 1,035 cases had genome sequences of sufficient quality to provide phylogenetic lineages. These genomes belonged to 26 distinct global lineages, indicating that there were multiple separate introductions into the region. Furthermore, 100 genetically-distinct UK lineages were detected demonstrating local evolution, at a rate of ~2 SNPs per month, and multiple co-occurring lineages as the pandemic progressed. Our analysis: identified a sublineage associated with 6 care facilities; found no evidence of reinfection in longitudinal samples; ruled out a nosocomial outbreak; identified 16 lineages in key workers which were not in patients indicating infection control measures were effective; found the D614G spike protein mutation which is linked to increased transmissibility dominates the samples and rapidly confirmed relatedness of cases in an outbreak at a food processing facility. The large-scale genome sequencing of SARS-CoV-2-positive samples has provided valuable additional data for public health epidemiology in the Norfolk region, and will continue to help identify and untangle hidden transmission chains as the pandemic evolves.